Stability and Symmetry Breaking in Metal Nanowires 
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A general linear stability analysis of simple metal nanowires is presented using a continuum ap- 
proach which correctly accounts for material-specific surface properties and electronic quantum-size 
effects. The competition between surface tension and electron-shell effects leads to a complex land- 
scape of stable structures as a function of diameter, cross section, and temperature. By considering 
arbitrary symmetry-breaking deformations, it is shown that the cylinder is the only generically 
stable structure. Nevertheless, a plethora of structures with broken axial symmetry is found at 
low conductance values, including wires with quadrupolar, hexapolar and octupolar cross sections. 
These non-integrable shapes are compared to previous results on elliptical cross sections, and their 
material-dependent relative stability is discussed. 
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I. INTRODUCTION 

Free-standing metal nanowires, suspended from electri- 
cal contacts at their ends, have been fabricated by a num- 
ber of different techniques. 1 Metal wires down to a single 
atom thick were extruded using a scanning-tunneling mi- 
croscope tip. 2,3 Metal nanobridges were shown to "self- 
assemble" under electron-beam irradiation of thin metal 
films, 4 ' 5 ' 6 leading to nearly perfect cylinders down to four 
atoms in diameter, with lengths up to fifteen nanome- 
ters. Systematic studies of nanowire properties for a vari- 
ety of materials were carried out using the mechanically- 
controllable break junction technique. 7 ' 8 ' 9 ' 10 ' 11,12 As the 
ultimate nanoscale conductors, metal nanowires are of 
great interest for nanotechnology. 

A remarkable feature of metal nanowires is that they 
are stable at all. Most atoms in such a thin wire are at the 
surface, with small coordination numbers, so that surface 
effects play a key role in their energetics. Indeed, macro- 
scopic arguments comparing the surface- induced stress to 
the yield strength indicate a minimum radius for solid- 
ity of order ten nanometers. 13 Below this critical radius, 
plastic flow would lead to a Rayleigh instability, 14 absent 
some other stabilizing mechanism. 

A series of experiments on alkali metal 
nanocontacts 7 ' 8 ' 9 identified electron-shell effects as 
another key mechanism influencing nanowire stability. 
Energetically-favorable structures were revealed as peaks 
in conductance histograms, periodic in the nanowire ra- 
dius, analogous to the electron-shell structure previously 
observed in metal clusters. 15 A supershell structure was 
also observed, 8 in the form of a periodic modulation 
of the peak heights. Recently, electron-shell structure 
has also been observed for the noble metals gold, 10 ' 11 
copper, 12 and silver. 12 

A theoretical analysis using the nanoscale free-electron 
model 16 ' 17 (NFEM) found that nanowire stability is 
determined by the competition of these two key fac- 
tors, surface tension and electron-shell effects. Both 
linear 13 ' 18 ' 19 and nonlinear 20,21 stability of axially sym- 



metric nanowires were investigated. It was found that 
the surface-tension driven instability can be completely 
suppressed in the vicinity of certain "magic radii." 

It is well known in the physics of crystals and molecules 
that a Jahn- Teller deformation breaking the symmetry of 
the system can be energetically favorable. In metal clus- 
ters, Jahn- Teller deformations are also very common, 22 ' 23 
and most of the observed structures show a broken spher- 
ical symmetry. By analogy, it is natural to assume that 
for nanowires, too, a breaking of axial symmetry can be 
energetically favorable, and lead to more stable deformed 
geometries. 

Recently, as a first step towards a stability analysis of 
wires with a more general cross section, elliptical wires 
were examined within the NFEM. 24 The sequence of sta- 
ble cylindrical and elliptical nanowires allows for a con- 
sistent interpretation of experimental conductance his- 
tograms for alkali metals, including both the electronic 
shell and supershell structures. 25 

Note that while the experimental manifestations of 
electron-shell structure are similar in metal clusters and 
nanowires, the Jahn- Teller effect plays out quite differ- 
ently in these two systems. The fundamental difference is 
that surface effects tend to stabilize clusters, while they 
are the source of the Rayleigh instability in nanowires. 18 
Therefore, Jahn- Teller deformations of clusters are very 
common and typically rather small, 22,23 while they only 
occur for a minority of stable nanowires and can be rather 
large. 

Cylindrical and elliptical nanowires are special, in the 
sense that the classical electron dynamics in these struc- 
tures is integrable. An open question is whether integra- 
bility per se plays a special role in nanowire stability. On 
the one hand, the shell effect is enhanced in integrable 
systems, 26 which argues in favor of integrability. On the 
other hand, the Jahn- Teller effect is driven by the lifting 
of degeneracy due to symmetry breaking, and the degen- 
eracy of states with angular momentum ±fj,h about the 
axis of symmetry is broken to leading order by a pertur- 
bation of the radius 5r{<fi) oc cos(2/i<^>), where <f> is the az- 
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imuthal angle, which renders the dynamics chaotic. Fur- 
thermore, for the case fx = 1 (quadrupolar cross section), 
the surface-energy cost of the perturbation is somewhat 
less than that of an elliptical deformation with the same 
aspect ratio. 

In this article, we perform a general linear stability 
analysis for straight metal nanowires, including arbitrar- 
ily shaped, non-integrable cross sections. By consider- 
ing all symmetry-breaking deformations, it is rigorously 
shown that the cylinder is the only generically stable 
structure. We derive the complete stability diagram for 
cylinders, and discuss why axial symmetry is special — in 
the sense that 75% of the experimentally observed alkali 
metal nanowires are indeed cylindrical, and that propor- 
tion goes up with increasing conductance. 

Nevertheless, a wide range of structures with broken 
axial symmetry is found at low conductance values, most 
of which show a reduced relative stability compared to 
the axially symmetric ones. We examine wires with 
quadrupolar, hexapolar and octupolar cross sections, and 
compare these non-integrable shapes to the previous re- 
sults on elliptical cross sections. Their relative stability 
and their degree of deformation are material-dependent, 
and we examine these properties for the alkali and noble 
metals. 

This paper is organized as follows: The assumptions 
and features of the nanoscale free-electron model are 
summarized in Sec. II, followed by the presentation of 
the linear stability analysis in Sec. III. Cylindrical wires 
are examined in Sec. IV, whereas Sec. V considers non- 
integrable cross sections. The material dependence of our 
results is discussed in Sec. VI. Finally, we summarize and 
discuss our results in Sec. VII. The Appendix provides 
additional technical details of our calculations. 



II. NANOSCALE FREE ELECTRON MODEL 

Guided by the importance of conduction electrons in 
the cohesion of metals, and by the success of the jcllium 
model in describing metal clusters, 15,27 the nanoscale 
free-electron model 16,17 replaces the metal ions by a uni- 
form, positively charged background that provides a con- 
fining potential for the electrons. The electron motion is 
free along the wire, and confined in the transverse di- 
rections. Due to the excellent screening 28,29 in metal 
wires with G > Go, where Go = 2e 2 /h is the con- 
ductance quantum, electron-electron interactions can in 
most cases be neglected. The surface properties of vari- 
ous metals can be fit by using appropriate surface bound- 
ary conditions. 24,30 

The NFEM is especially suitable for alkali metals, 
but is also adequate to describe shell effects due to the 
conduction-band s-electrons in other monovalent metals, 
such as noble metals, and in particular, gold. The ex- 
perimental observation of a crossover from atomic-shell 
to electron-shell effects with decreasing radius in both 
metal clusters 31 and nanowires 9 justifies a posteriori the 



use of the NFEM in the later regime. 

Since a nanowire connecting two macroscopic elec- 
trodes is an open quantum system, the Schrodingcr equa- 
tion is most naturally formulated as a scattering problem. 
The fundamental theoretical quantity is the scattering 
matrix S(E) connecting incoming and outgoing asymp- 
totic states of conduction electrons in the electrodes. 
Thermodynamic properties can be expressed in terms of 
the scattering matrix through the electronic density of 
states 32 

from which the relevant thermodynamic potential for an 
open system, namely the grand canonical potential ft, is 
obtained as 



dE D(E) In 



1 + e 



-0{E-n) 



(2) 



Here (3 = (ksT) -1 is the inverse temperature and \i is the 
electron chemical potential, specified by the macroscopic 
electrodes. Eqs. (1) and (2) include a factor of 2 for spin 
degeneracy. 

Any extensive thermodynamic quantity can be ex- 
pressed as the sum of a Weyl expansion, which depends 
on geometrical quantities such as the system volume V, 
surface area 5, and integrated mean curvature C, and an 
oscillatory shell-correction due to quantum-size effects. 26 
In particular, the grand canonical potential (2) can be 
written as 



n = -ojV + ct s s -j s c + sn, 



(3) 



where the energy density —u>, surface tension coefficient 
<7 S , and curvature energy — 7 S arc, in general, material- 
and temperature-dependent coefficients. On the other 
hand, the shell correction Sfl can be shown, based on very 
general arguments, 29,33,34 to be a single-particle effect, 
which is well described by the NFEM. 

The ionic degrees of freedom in the NFEM are mod- 
eled as an incompressible fluid. 13,20,21 This takes into ac- 
count, to lowest order, the hard-core repulsion of the core 
electrons as well as the exchange energy of the conduc- 
tion electrons. Any atom-conserving deformation of the 
structure is therefore subject to a constraint of the form 



J\f = kpV - r] k 2 F S + £ k F C = const, 



(4) 



where kp is the Fermi wavevector. The parameters ij and 
£ can be adjusted so as to fix the values of the effective 
surface tension and curvature energy in the free-electron 
model to the material-specific values, as discussed in de- 
tail in Sec. VI. However, the shell correction 6£l is (to 
leading order) independent 24,34 of r\ and £. 

In this article, the shell correction SQ is calculated 
quantum-mcchanically, from Eqs. (1) and (2). Nonethe- 
less, in interpreting our quantum-mechanical results, it is 
useful to appeal to semiclassical arguments, 26,35 accord- 
ing to which the shell correction can be computed by a 
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Gutzwillcr-type sum 36 over classical periodic orbits, that 
takes the form of an integral over the length L of the 
nanowire 

E 

Sn = -^ dzS2A po (z,T)cos(k F L po (z) + 9 po ), (5) 
Xf J o po 

where Ef and Xf are the Fermi energy and Fermi wave- 
length, respectively. Here the sum includes all classi- 
cal periodic orbits lying in the plane of the wire's cross 
section, L po is the orbit length, and 9 po is an orbit- 
dependent phase. Semiclassically, 26 the temperature de- 
pendence of 6Q appears only in the dimcnsionlcss ampli- 
tude A po , which also depends on the stability, symmetry, 
and length of the periodic orbit. The orbit amplitude has 
the form 26 

where T F = E F /k B is the Fermi temperature, and 
v po = 1 for a continuous family of degenerate orbits, 
which arise when the cross section has the form of an intc- 
grable billiard, such as a disk or an ellipse, and v po = 3/2 
for an isolated orbit in a wire with a cross section whose 
dynamics is fully chaotic (see Ref. 37 for details). 

It is useful to define the r.m.s radius p of the cross sec- 
tion, which is related to the cross-sectional area by A = 
up 2 . Then L po oc p, so that <5£1 <~ p^ 1 for an integrable 
cross section, and decays exponentially for temperatures 
exceeding the radius-dependent characteristic tempera- 
ture T p = Tp/kFP- As shown below, these semiclassical 
estimates of the temperature- and size-dependence of <5f2 
are in accord with our quantum- mechanical results. 

III. LINEAR STABILITY ANALYSIS 

Let us consider a wire with uniform cross section 
aligned along the z-axis. Its surface is given by the ra- 
dius function r = R(<p; p, A) in cylindrical coordinates 
r, ip, z, where the parameters {p, A} describe the cross 
section geometry, chosen such that the cross-sectional 
area is given by A = Trp 2 , whereas the shape is deter- 
mined by a set of dimensionless parameters, composing 
the vector A. Without loss of generality, we can associate 
A = with a cylindrical wire. 

A small z-dependent perturbation of a wire of length 
L and initial cross section (p, A) can be written in terms 
of a Fourier series as 

p(z) = p + eSp(z) = p + e^p q e lq \ 

i 

\(z) = A + eSA{z) = A + e^A q e iqz , (7) 

i 

where the dimensionless small parameter e sets the size 
of the perturbation. Assuming periodic boundary con- 
ditions, the perturbation wave vectors q must be integer 



multiples of 2w/L. In order to ensure that p(z) and A(z) 
are real, we have p- q = p* and A_ g = A*. 

The structural stability of metal nanowires is 
governed by the response to long-wavelength 
perturbations, 13 ' 17,18 ' 24 while the response to short- 
wavelength perturbations controls surface quantum 
fluctuations in long wires. 19,38 Since we consider free 
electrons, the stationary Schrodinger equation is given 
by — = (2m e E /h 2 )^ , where A is the Laplace oper- 
ator in cylindrical coordinates. In the long-wavelength 
limit, the Schrodinger equation may be solved in the 
adiabatic approximation, for which transverse and longi- 
tudinal motions decouple. Therefore, we use the ansatz 
^(r, <p, z) — x{ r , Vi z )^{ z )i an d neglect all ^-derivatives 
of the transverse wavefunction x, so that the Schrodinger 
equation becomes 

f d 2 Id Id 2 2m e _ , A . , . , 

{^ + r^W + ^ En{Z) ) Xn{r ^ Z) = °' (8) 
/ f> 2 2m \ 
\dz^ + -ff(E-E n (z)))$ n (z) = 0.(9) 

First, the transverse problem (8) is solved at fixed z in 
order to determine the transverse eigenenergies E n (z). 
With our choice of cross section parametrization, their 
dependence on geometry can be written as 

where the function 7«(A) remains to be determined. For 
an axisymmctric wire, the 7„ are given by the roots of the 
Bessel functions, whereas for an elliptical cross section, 
the 7„ are given by the roots of the modified Mathieu 
functions. 24 In general, the shape-dependent 7«(A) are 
determined numerically (sec Appendix for details). 

We are then left with a series of effective one- 
dimensional scattering problems (Eq. 9) for the longi- 
tudinal wave functions <f> n (z), in which the transverse 
eigenenergies E n i y p{z), A(z)) act as additional potentials 
for the motion along the wire. These scattering prob- 
lems can be solved using the WKB approximation 16 : The 
grand canonical potential is given by 

n[T;p,A} = J™ dE(-^)E[E;p,A], (11) 

where / = (l+exp[(£— ^)/fcgT]) _1 is the Fermi function 
at temperature T and chemical potential p w E F , and S 
is given by 

~, = -^ ? (^i^) s % 2 , 

The sum runs over all open channels n, for which 
E n (p,A)<E. 

The energetic cost of a small deformation of the wire 
can be calculated by expanding Eq. (11) as a series in 
the parameter e, 

n = fi<°> +ef!« +e 2 n^ +0(e 3 ). (13) 
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A nanowire with initial cross section (p, A) is energeti- 
cally stable at temperature T if and only if fiW (p, A) = 
and Q( 2 )(p,A) > for every possible deformation 
(Sp,SA) satisfying the constraint (4). 

A straightforward expansion of Eq. (1 1) at T = yields 



L/X f 



E„ 



Ao-Va-Ek — 2E n -^- 
P 



= ^E 



Pq/p\ f A pp A pA \ f Pq/P 



A q J \-A-Ap Aaa J \ A 



(14) 



(15) 



where the elements of the matrix A in Eq. (15) are given 
by 



Ap P — 



E 



AE„ 



E 



3/2 



3y/E F -E n 



E 



4£: 



E 



3/2 



= E 



E 



3/2 



VE F -E n - 
2E'^E F -E n 



E n 



2sjE F ~E n _ 

K ■ {E'J 



,(16) 



E n {p,A), E' n = V A E n 
with respect to A, 
is the matrix of second derivatives, all 



\J Ef — E n 

p A denotes the 
and = 



Here E n 
gradient of E, 
d 2 E n /dA i dA j \.- A 

evaluated at (p, A). 

The number of independent Fourier coefficients in Eq. 
(7) is restricted through the constraint (4) on allowed 
deformations. Hence, after evaluating the change of the 
geometric quantities V, S, and C due to the deformation, 
we can use Eq. (4) to express po in terms of the other 
Fourier coefficients. This yields an expansion p a = p^ + 
C(e 2 ),with 



e Pb 

(o) 
Ph 



>1 



k F p — f}P 
1 

2{k F p-r 1 P) 



P'-A 



o ■ 



(17) 



E 



A q ) {v(P') T 



T)P' 
X)P" 



Pq/P 
A, 



where the function P = P(A) is related to the perimeter 
V by V(p, A) = 2npP, the first derivative P' = V A P| A is 
a vector, and the second derivative P" = d 2 P/dAidAj\ A 
is a matrix. Note that in the long-wavelength limit con- 
sidered here, Eq. (17) is independent of the parameter £ 
appearing in the constraint (4), so that our results are 
independent of the value of £. Inserting the expansion 
for p in Eq. (14), the first-order change of the energy 
under a constant Af perturbation is given by 



A 



— 

L 



AT 



fEl 



VEf 



E n 



(18) 



EL 



2r]E n 



k,Fp — TfP 



P' ■ A , 



and the second-order term is given by Eq. (15), with the 
matrix A replaced by 



lpp " App+ k F p r -vP^ 



Ak F p E n \/ Ep — E n 




A\\ = A AA 



kFP — TfP 



The stability condition, ^- 2 \p, A)\^ > 0, requires that 
the stability matrix A be positive definite. The results at 
finite temperature are obtained in a similar fashion, by 
integrating Eq. (11) numerically. 

IV. GENERAL STABILITY OF CYLINDERS 

As a first application of the method presented in the 
previous sections, we derive the complete stability dia- 
gram for cylinders, i.e., we determine the radii of cylindri- 
cal wires that are linearly stable with respect to arbitrary 
small, long-wavelength deformations. We choose the fol- 
lowing general form of the radius function describing the 
surface of the wire, 

\ V rn m / 

(20) 

The sum runs over a set IM of positive integer values, 
and the deformation parameters A m governing deviations 
from axial symmetry are such that J2 m < 2. Note 
that the dipolc deformation (m = 1) corresponds, in 
leading order, to a simple translation, plus higher-order 
multipole deformations. Therefore we can restrict our 
analysis to to > 1. 

At first sight, considering arbitrary deformations, and 
therefore theoretically an infinite number of perturbation 
parameters {\ m ,(p m } seems a formidable task. Fortu- 
nately, we find that the stability matrix A for a cylinder 
is diagonal, and therefore the different Fourier contribu- 
tions of the deformation decouple. 

This can be seen within perturbation theory for small 
deviations from axial symmetry, derived in the Appendix: 
for energy levels n that are non-degenerate at A = 0, i.e. 
levels with orbital quantum number p n = 0, we obtain 



asn(o) 

dA 



and 



d 2 E n (0) 
dAidAj y ' 



(21) 



whereas for two energy levels n and n' that are degenerate 
at A = 0, i.e. those for which p n = —p n > ^ 0, we get 



dE n (0) dE n ,(0) 



dA 



dA 



(22) 
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FIG. 1: Stability of metal 
nanocylinders versus electrical 
conductance and temperature. 
Dark gray areas indicate sta- 
bility with respect to arbitrary 
small deformations. Tempera- 
ture is displayed in units of T p = 
Tf/UfP (see text). The surface 
tension was taken as 0.22 N/m, 
corresponding to Na. 39 For com- 
parison, the upper panel also 
shows wires that are stable to- 
wards axisymmetric deforma- 
tions only (light gray areas), but 
which are unstable when allow- 
ing symmetry breaking defor- 
mations. The lower panel in- 
cludes a plot of the activation 
energy AE (solid line), repro- 
duced from Ref. 21, which de- 
termines the nanowire lifetime 
T = r exp(A£/fc s T). 



and 



d 2 E n (0) d 2 E n ,(0) 
— H — oc da 

dAidAj dkidkj 3 



(23) 



Since the matrix A, given by Eq. (16), includes a sum 
over all open channels, the relations (21-23) imply that 
A is diagonal. Furthermore, one can straightforwardly 
show, by expanding the perimeter function P in a se- 
ries in A, that P' = Va-P|a=o = and the matrix 
P" = d 2 P/dAidAj\A=o is diagonal as well. Therefore, 
we conclude that the stability matrix A, including the 
corrections from the constraint (4), is diagonal. This 
also shows that the first order correction (18) to fi is 
identically zero for cylinders. This result implies a great 
simplification of the problem, since it allows to deter- 
mine the stability of cylindrical wires with respect to ar- 
bitrary deformations through the study of a set of pure 
m-deformations, i.e. deformations as given by Eq. (20) 
with only one non-zero X m . 

Fig. 1 shows the stable cylindrical wires in dark gray 
as a function of temperature. The surface tension was 
fixed at the value 39 for Na, a s = 0.22 N/m (i.e. r\ = 1.1). 
The temperature is given in units of the characteristic 
temperature T p = Tf/UfP, reflecting the temperature 
dependence of the shell correction (5) to the wire energy: 
Since stability is determined by a competition between 
the surface and shell contributions, 17,18 the maximum 
temperature of linear stability scales with T p , the sur- 
face contribution being only weakly dependent on tem- 
perature. The £-axis is given by the corrected Sharvin 
conductance 



Gs = G 



k 2 F p 2 k F P 

4 47T 



(24) 



where V is the perimeter of the cross section. 

The stability diagram was obtained by intersecting a 
set of individual stability diagrams allowing cos(m</?)- 
deformations with m = 2, ...,8. Multipolar deformations 
become increasingly costly with increasing m, their sur- 
face energy scaling as pm 2 , and we find that including 
m > 6 deformations does not further modify the final 
result. 

Figure 1 confirms the extraordinary stability of a set 
of wires with so called "magic radii," which were identi- 
fied in previous studies of stability under axisymmetric 
perturbations alone. 13,18 ' 19,21 They exhibit conductance 
values G/G = 1, 3, 6, 12, 17, 23, 34, 42, 51,... How- 
ever, some wires that are barely stable when considering 
only axisymmetric perturbations, 13,18 ' 19 e.g., G/Go = 5, 
10, 14,..., shown in light gray in the top panel of Fig. 
1, are found to be unstable when allowing more general, 
symmetry-breaking deformations. 

The heights of the dominant stability peaks in Fig. 
1 exhibit a periodic modulation, with minima occurring 
near G s /G = 9, 29, 59, 117, 170, and 255. The posi- 
tions of the first four minima are in perfect agreement 
with the observed supershell structure in conductance 
histograms of alkali metal nanowires, 8 while we predict 
the next two minima of the scries. The supershell struc- 
ture is expected to continue at larger conductance, but 
may be obscured by atomic-shell effects in experiments. 
Interestingly, the nodes of the supershell structure, where 
the shell effect for a cylinder is suppressed, are precisely 
where the most stable deformed nanowires are predicted 
to occur (see Ref. 24 and the discussion below). Thus 
symmetry breaking distortions and the supershell effect 
are inextricably linked. 

Linear stability is a necessary — but not a sufficient — 
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(a) 



(b) 



K» 0.4 




FIG. 2: (color online) (a) 
Linearly stable quadrupo- 
lar (A2 > 0) Na wires 
[thick lines] at tempera- 
ture T = 0.03 T p , shown in 
the configuration space of 
the two parameters A2 and 
p describing the quadrupo- 
lar shape, R(f) = 
p(V'l-A2/2 + A 2 cos(2v9)). 
p is related to the 
Sharvin conductance 
by Eq. (24), which gives 
y/Gs/Go ~ k F p/2. The 
thin gray lines show the 
thresholds for the opening 
of new channels. Some 
geometries are shown on 
the right axis, illustrating 
their deviations from axial 
symmetry; (b) Stability 
diagram for quadrupolar 
Na wires. The gray(color) 
scale reflects the value of 
A2 for the stable wires. 



condition for a nanostructure to be observed experimen- 
tally. The linearly stable nanocylinders revealed in the 
above analysis are in fact metastable structures, and 
an analysis of their lifetime has recently been carried 
out within an axisymmetric stochastic field theory. 21 
In the lower panel of Fig. 1, the zones of linear sta- 
bility are superimposed on the activation energy AE 
(solid line), calculated in Ref. 21, which determines the 
nanowire lifetime r through the Kramers formula r = 
To exp(A-E//csT). As can be seen, there is a strong cor- 
relation between the size of the activation barriers and 
the height of the stable fingers in the linear stability anal- 
ysis. This suggests that the linear stability analysis, with 
temperature expressed in units of T p = Tp/kFp, provides 
a good measure of the total stability of metal nanowires. 
In particular, the "universal" stability 21 of the most sta- 
ble cylinders is reproduced, wherein the absolute stability 
of the magic cylinders is essentially independent of radius 
(aside from the small supershell oscillations). 

For a given material, the typical value AE for the 
magic cylinders is universal, i.e. it depends only on the 
surface tension 21 cr, 



AE « 0.6\Jha s /m 



(25) 



where m e is the electron effective mass, which is of order 
the free electron mass in simple metals. When comparing 
the absolute stability for different materials, it is not suf- 
ficient to simply compare the heights of the stable fingers 
in Fig. 1, which depend only logarithmically on the ratio 



a s /Epkp through Eqs. (3) and (5); the a^J 2 dependence 
of the activation barrier, which determines the lifetime, 
must be taken into account explicitly (see Sec. VI). 

Since the activation barriers are calculated asymptot- 
ically for low temperatures, 21 they exhibit all of the fine 
structure of the low-temperature shell correction. This 
fine structure is also seen in the linear stability analysis 
at low temperature (see Fig. 1, lower panel). However, 
at higher temperatures, the fine structure in the shell po- 
tential is smoothed out, as described in Eqs. (5) and (6), 
and some of the separate low-temperature stability zones 
merge into arches at higher temperatures. Prominent ex- 
amples of this phenomenon, the unusual thermodynamics 
of which is discussed in Ref. 13, occur at G5/G0 ~ 67, 
80, and 255. 



V. BREAKING AXIAL SYMMETRY 

In a recent paper, 24 the stability of wires with elliptic 
cross sections was determined. Elliptic deformations are 
special, in that the system remains integrable. Here we 
focus on wires with a cos(m(p)-deformed cross section (i.e. 
having m-fold symmetry), a special case of Eq. (20) with 
only one non-zero A m , which renders the system non- 
integrable. In the following, we will discuss quadrupolar 
deformations with m — 2, which are the most energeti- 
cally favorable of the multipole deformations, as well as 
the higher-order multipoles with m < 6. Deformations 
with higher m cost significantly more surface energy, and 
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thus yield fewer stable configurations. The stability anal- 
ysis of wires with 5- and 6-fold symmetry reveals only a 
few barely stable configurations (besides the axisymmet- 
ric wires discussed above), whereas for to > 6 no sta- 
ble deformed wires were found at all in the temperature 
range considered. 

Figure 2(a) shows the stable configurations [thick 
curves] for quadrupolar Na wires at temperature T = 
0.03 T p . These are determined by the intersection of the 
stationary curves, Vt^(p, \2)\n — 0, and the convex re- 
gions, fi( 2 )(p, A 2 )|x > 0. Also shown are the thresh- 
olds to open new conducting channels, i.e., the energy 
eigenvalues of the two-dimensional billiard comprising 
the cross section of the wire. Note that the stable struc- 
tures lie in the gaps of the spectrum of channel thresh- 
olds, as is the case for axisymmetric wires. 18 

Combining results for all temperatures yields the sta- 
bility diagram shown in Fig. 2(b). Note that the max- 
imum deformation of the stable structures decreases 
strongly with increasing conductance. Nanowires with 
highly- deformed cross sections are only stable at small 
conductance, unlike the magic cylinders, which are pre- 
dicted to occur for arbitrarily large conductance (until 
the transition to crystalline structures at sufficiently large 
radius). 

The most stable wires are found at the same conduc- 
tance values, and with similar semi-axis ratios, as the 
stable elliptical wires, discussed in Ref. 24. This is not 
surprising, since the ellipse can be approximated to lead- 
ing order by a quadrupole, 

Re\i{<p) ~ (1 - k 2 ) + (2k + 7k 3 ) cos(2^) + 3k 2 cos(4^) 
+5k 3 cos(6^) +£>(k 4 ), (26) 

where k = |(e 2 — l)/(e 2 + 1) is given in terms of the 
aspect ratio e. The maximum temperature T max up to 
which the wires are stable is in general of the same order 
as found for elliptical wires. Nevertheless, for the stable 
wire at G = 2Go, T max is 20 percent larger compared to 
the corresponding elliptical wire, indicating that a wire 
with two conducting channels will show a peanut-shaped 
cross section, rather than an elliptic one. Remarkably, 
this result, which was obtained by minimizing the elec- 
tronic energy — and thus contains nothing of the atomic 
structure — is exactly what one would expect for a wire 
with two atoms in the cross section. 

In order to shed light on the possible competi- 
tion between integrable (elliptical) and non-integrable 
(quadrupolar) shapes, we analyzed the stability of wires 
with simultaneous to = 2 and m — 4 deformations, which 
can interpolate between elliptical and quadrupolar cross 
sections [cf. Eqs. (20) and (26)]. The results are summa- 
rized in Table I. For large deformations, the quadrupole 
is more stable, due to its lower surface energy, but no 
overall preference for a special geometry was found; sta- 
ble wires with both positive and negative values of A4 
exist. The values of A4 were found to be small, indicat- 
ing that the addition of further deformations has little 
influence on the stability diagram. Remarkably, the wire 



(a) G. / a s 




0.05 0.1 0.15 



FIG. 3: (color online) Linear stability of cos(m(/p)-deformed 
wires for (a) m = 3, (b) m = 4, (c) m — 5, and (d) m = 6, 
as a function of temperature (in units of T p ). A logarithmic 
scale is used for the temperature since wires with m > 4 are 
much less stable than those for m < 3. A sketch of the cross 
section is shown on the right. The deformation parameter 
A m is encoded via a gray(color) scale, common to all four 
diagrams. In the interest of clarity, only non-cylindrical wires 
(with A m > 0) are shown. The surface tension was chosen to 
represent Na. 

with G = 2Gq is even closer to two touching cylinders 
when 4-fold deformations are added to the quadrupolar 
ones. 

The linearly stable wires with 3-, 4-, 5-, and 6-fold 
symmetry are shown in Fig. 3 as a function of tempera- 
ture. For clarity, the stable cylindrical wires are omitted, 
and a logarithmic scale is used for the temperature axis. 
The value of the deformation parameter A m is encoded 
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G 






Ellipse 




A2 cos(2ip) 


A2 cos(2</j) + A4 cos(4</s) 


[Go] 


£ 


A 2 


A 4 




e 


A 2 




A 2 


A 4 




2 


1.65 


0.24 


0.0440 


0.40 


1.7 


0.26 


0.50 


0.25 


-0.04 


0.50 


5 


1.32 


0.14 


0.0140 


0.44 


1.32 


0.14 


0.49 


0.15 


-0.03/+0.08 


0.49/0.23 


9 


1.24 


0.11 


0.0086 


0.50 


1.22 


0.10 


0.50 


0.10 


0.03/-0.03 


0.48/0.44 


29 


1.14 


0.07 


0.0032 


0.54 


1.13 


0.06 


0.54 


0.07/0.05 


0.008/-0.015 


0.44/0.44 


59 


1.09 


0.04 


0.0014 


0.50 


1.11 


0.05 


0.49 


0.04 


-0.01 


0.49 


72 


1.08 


0.04 


0.0011 


0.40 


1.08 


0.04 


0.39 


0.04 


0.009 


0.36 



TABLE I: Comparison of the most stable deformed wires with elliptic (columns 2-5), quadrupolar (columns 6-8), and more 
general (columns 9-11) cross sections. The first column gives the quantized conductance of the corresponding wire. For the 
elliptic and quadrupolar wires, both the aspect ratio and the value of the deformation parameters (obtained from Eq. (26) for 
the former) are given. The maximum temperature of stability T max is given for each wire. The values for the elliptic wires 
are taken from Ref. 24. Note that multiple values indicate multiple stable wires with the same conductance. In all cases the 
surface tension was set to 0.22 N/m, corresponding to Na. 39 



through a gray(color) scale, which is common to all sta- 
bility diagrams shown in this figure, but is different than 
the one of Fig. 2. Compared to the quadrupolar wires, 
the number of stable configurations, the maximum tem- 
perature of stability, and the size of the deformations 
involved, all decrease rapidly with increasing order m of 
the deformation. For m > 6, we no longer find stable 
geometries in the temperature range considered. All this 
reflects the increase in surface energy with increasing or- 
der m of the deformation. 

From Fig. 3(b), we can extract a series of stable wires 
with four-fold symmetry, for which the maximum defor- 
mation of A4 ~ 0.1 is found at a conductance of G = 4 Go- 
Although this geometry is expected to be far less stable 
than the neighboring cylindrical and quadrupolar wires, 
it has nevertheless likely been observed as a shoulder in 
Na conductance histograms. 25 

VI. MATERIAL DEPENDENCE 

Within the NFEM there is only one parameter enter- 
ing the calculation apart from the contact geometry: the 
Fermi energy Ep, which is material dependent and in 
general well known (see Tab. II). Nevertheless, we have 
seen in the previous section that the energy cost of a de- 
formation due to surface and curvature energy, which can 
vary significantly for different materials, plays a crucial 
role in determining the stability of a nanowire. Using the 
NFEM a priori implies the macroscopic free energy den- 
sity oj = 2Epkp/15n 2 , the macroscopic surface energy 
a s = Ephp/lfm, and the macroscopic curvature energy 
7 S = 2Epkp/9ir 2 . When drawing conclusions for met- 
als having surface tensions and curvature energies that 
are rather different from these values, we have to think 
of an appropriate way to include these material-specific 
properties in our calculation. 

Other authors have calculated the surface energy us- 
ing a free-electron model with a rectangular confining 
potential of finite height, 40 ' 41 depending on the electron 
work function, without being able to reproduce the values 



found in experiments. The agreement does not improve 
even when using a self-consistent confining potential. 41 
When working with a free-electron model, contributions 
of correlation and exchange energy are not included, 
which are found to play an essential role for a correct 
treatment of the surface energy. A discussion of the cal- 
culation of the surface energy of a jellium metal beyond 
the free-electron model can be found in Ref. 42. 

A convenient way of modeling the material properties 
without losing the pleasant features of the NFEM is via 
the implementation of the constraint (4) on the defor- 
mation, which interpolates between volume conservation 
(?7 = £ = 0) and treating the semiclassical expectation 
value for the charge N Weyl as an invariant (77 = 37r/8, 

£=!)• 

Consider the grand canonical potential of a free- 
electron gas confined within a given geometry by hard- 
wall boundaries, as given by Eq. (3). The change in en- 
ergy due to a deformation is 

Afl = -cuAV + a s AS- 7s AC + A[5Sl] 

-( 7i - ^)AC + A[«l], 

where we have used the constraint (4) to eliminate V. 
Now the prefactors of the change in surface AS and the 
change in integrated mean curvature AC can be iden- 
tified as effective surface tension and curvature energy, 
respectively. They can be adjusted to fit the material 
properties by appropriate choice 43 of the parameters r? 
and £ (see Tab. II). As noted previously, results in the 
long- wavelength limit are independent of £. 

The stability diagrams discussed in the previous sec- 
tions all represent sodium nanowires. Results for other 
s-orbital metals are similar in respect to the number of 
stable configurations and the conductance of the wires. 
On the other hand, the deviations from axial symmetry 
and the relative stability of Jahn- Teller deformed wires 
is sensitive to the material specific surface tension a s 
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Element 


Li 


Na 


K 


Cu 


Ag 


Au 


E F [eV] 
kF [rim -1 ] 


4.74 
11.2 


3.24 
9.2 


2.12 
7.5 


7.00 
13.6 


5.49 
12.0 


5.53 
12.1 


a a [meV/A 2 ] 
o a [E F k 2 F ] 
V 


27.2 
0.0046 
1.135 


13.6 
0.0050 
1.105 


7.58 
0.0064 
1.001 


93.3 
0.0072 
0.939 


64.9 
0.0082 
0.866 


78.5 
0.0097 
0.755 


7 S [meV/A] 
7s [E F k F ] 
t 


62.0 
0.0117 
0.802 


24.6 
0.0082 
1.06 


14.9 
0.0094 
0.971 


119 
0.0125 
0.741 


96.4 
0.0146 
0.583 


161 
0.0240 
-0.111 



TABLE II: Material parameters of several monovalent metals: 
Fermi energy Ef, 44 Fermi wavevector kF, 44 surface tension 
o s , and curvature energy 7s, 39 along with the corresponding 
values of r\ and £. 



and Fermi temperature Tp. The relative stability of the 
highly deformed wires decreases with increasing surface 
tension a s /E F k F , measured in intrinsic units, and this 
decrease becomes stronger with increasing order m of the 
deformation. Therefore, for the simple metals under con- 
sideration (Tab. II), deformed Li wires have the highest, 
and Au wires have the lowest relative stability compared 
to cylinders of "magic radii." 

However, concerning absolute stability, one must con- 
sider that the lifetime of a metastable nanowire also de- 
pends on the surface tension [cf. Eq. (25)]. Thus, al- 
though the deformed wires have reduced relative stabil- 
ity in the noble metals, this effect is compensated by the 
greater absolute stability of these materials. 

In Fig. 4, the linearly stable quadrupolar wires for gold 
are shown. The vertical scale was adjusted according 
to the arguments discussed above, in order to allow a 
comparison with the absolute stability for Na wires, as 
shown in Fig. 2. 



VII. SUMMARY AND DISCUSSION 

In this article, we have performed a complete linear 
stability analysis of metal nanowires within the NFEM. 
That is to say, we have examined the stability with re- 




FIG. 4: Linearly stable Au wires with quadrupolar cross sec- 
tion. The deformation parameter A2 > is coded via a 
color(gray) scale. 



spect to small, long-wavelength perturbations of metal 
nanostructures of arbitrary cross section, with transla- 
tional symmetry along an axis perpendicular to the cross 
section. We have considered both cross sections for which 
the electron motion is integrable, namely circles and el- 
lipses, and nonintegrable cross sections, described by a 
multipole expansion (20). 

This general stability analysis confirms the central 
conclusions of previous analyses, which were limited to 
circular 13 ' 18,19 and elliptical cross sections. 24,25 First, the 
existence of a sequence of "magic" cylindrical wires of ex- 
ceptional stability, with electrical conductance G/Gq — 
1, 3, 6, 12, 17, 23, 34, 42, 51,..., was confirmed. Second, 
the existence of a series of very stable deformed wires, 
with two-fold symmetric cross sections, which occur at 
the nodes of the cylindrical shell structure, at G/Go = 2, 
5, 9, 29, 59,..., was verified. For the thinnest of these 
wires (G/Go = 2, 5), which have the largest deforma- 
tions, the quadrupolar cross section was found to be more 
stable than the ellipse, while the cross sections of the 
thicker wires in the sequence were found to be essentially 
elliptical. In addition, a number of new stable struc- 
tures with quadrupolar, hexapolar, and octupolar cross 
sections have been identified. 

Perhaps the most surprising result of our general linear 
stability analysis is the special role played by cylinders. 
Roughly 75% of the principal structures observed in con- 
ductance histograms 8 of alkali metals correspond to the 
magic cylinders. The remaining 25% are comprised of 
wires with elliptical/quadrupolar cross sections. Despite 
expanding the configuration space of allowed structures 
significantly in the present analysis, no new structures of 
comparable stability have been identified for the alkali or 
noble metals. The role of symmetry in the stability of 
metal nanowires is thus fundamentally different from the 
case of atomic nuclei 45 or metal clusters, 22,23 where the 
vast majority of stable structures have broken symmetry. 

The crucial difference between the stability of metal 
nanowires and metal clusters is not the shell effect, 
which is similar in both cases, but rather the surface en- 
ergy, which favors the sphere, but abhors the cylinder. 14 
In order to understand the unique stability of metal 
nanocylinders, in contrast to other wire structures, it is 
useful to consider the semiclassical size of the competing 
surface and shell contributions: From Eq. (3), the surface 
energy scales as 



a, 



0(k FP ), 



where p is the r.m.s. radius of the wire, 
hand, Eq. (5) implies that 



(28) 



On the other 



EpkpL 



o(k FP y 



(29) 



where v = 1 for circular or elliptical cross sections, and 
3/2 > v > 1 for nonintegrable cross sections. For a 
typical metal nanowire in the domain of validity of the 
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NFEM, k F p ~ 10 (i.e., G s /G ~ 20). Thus, the term ft s 
which gives rise to the Rayleigh instability ' 18 is roughly 
two orders of magnitude larger than the shell correction 
SO, which stabilizes the wire — even in the best-case sce- 
nario of an integrable cross section! (The volume contri- 
bution to 0, though larger still than O s , does not play a 
role in deformations that conserve the number of atoms.) 
How then are we to understand the stability of metal 
nanowires? 

Stability is not determined by the energy itself, but by 
the stationarity and convexity of the energy functional. 
Because <50 is a rapidly oscillating function [cf. Eq. (5)], 
its derivatives are large in a scmiclassical sense. 46 The 
first and second variations of SO scale as 



Epk F h 



o(k FP y 



so^ 

EpkpL 



0{k FP f 



(30) 



Thus, for the case of an integrable cross section (y = 1), 
SO. 1 - 2 ) is of the same semiclassical size as O^. However, 
SO^ is always semiclassically smaller than 0^\ The 
stationarity condition 0^ = is more difficult to satisfy 
than the convexity condition 0^ > 0, which can be 
satisfied for any integrable cross section. 

Cylinders are special because O^ = by symmetry, so 
stability is determined solely by convexity, where surface 
and shell effects compete on an equal footing (cf. Fig. 
1). In contrast, non-axisymmetric shapes can only occur 
for small kpp. More precisely, the maximum deviation 
from axisymmctry is a decreasing function of kpp (see 

Figs. 2 and 3). Both 50^ and O^ vanish as A — > 0, but 
whereas SO 1 - 1 ) is an oscillatory function of A [cf. Eq. (5)], 
O^ grows smoothly with increasing A: 

n^/iw a s \ m (m 2 - l)2wp. (31) 

The stationarity condition O^ = can therefore only 
be satisfied for 



A m < 



const. / E F k 2 F 



1 



(k F pY 



(32) 



Indeed, the decrease in the maximum deformations as a 
function of k F p in Figs. 2, 3, and 4 is consistent with Eq. 
(32) with 3/2 > v > 1. 

Thus cylinders are the only gcnerically stable struc- 
tures, and account for the principal peaks in conductance 
histograms of monovalent metals, while structures with 
broken symmetry are favored only at the nodes of the 
cylindrical supershell structure. 



APPENDIX: ENERGY SPECTRUM OF A 
2D-BILLIARD 

The essential input for the stability analysis of a 
nanowire (Sec. Ill) is the energy spectrum of a two- 
dimensional billiard for which the boundary (in polar co- 
ordinates r, ip) is given in terms of the radius function 
R(<p;p,A) = R v of the wire, by the condition r = R v . 
This provides the transverse eigenenergies E n [p, A) en- 
tering the expressions (18) and (19). It is helpful to 
introduce a change of coordinates, f — r/R v , which 
simplifies the boundary to a circle of unit radius. The 
Hamiltonian in the new coordinates (f, ip) can be writ- 
ten as H = JT(°) + SH, where = -R~ 2 A f ^ is the 
Hamiltonian for an axisymmetric problem and Af j¥ , is 
the Laplace operator. For simplicity of notation we use 
units in which h 2 /2m e = 1. On the other hand, SH con- 
tains the terms due to deviations from axial symmetry 
and is found to read 



SH = - 



K o 2 2 i?;i» 2% 

df 2 i?3 f df dip 



f df' 



(A.l) 



where a prime indicates differentiation with respect to 
ip. The energy spectrum can be obtained by numerical 
diagonalization of H in an appropriate basis, that we 
choose to be 



1 



e^Ml^f), (A.2) 



where the are the roots of the Bessel functions J^. 
We simplify the notation by using roman indices j as 
multi-indices representing a set of quantum numbers 
(fij,Vj). Note that the states (A.2), though normalized, 
are not orthogonal, since the scalar product in the new 
coordinates, 



U\k) 



/ R%d<p 
Jo Jo 



fdr**Af,<p)* k {r,<p), (A.3) 



involves a factor R 2 . However, the {^j} still form a ba- 
sis and the energy spectrum can be obtained by solving 
a generalized eigenvalue problem. Defining the matri- 
ces Hjk = (j | H\ k) and Bjk = (j\k) we have to solve 
det[W — EB] = 0. Taking matrix elements of SH we ob- 
tain after some algebra, 
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(j\6H\k) = (j 



j ''-j - 1 



+ (Pj+Pk)(j 



2 2 



vhl i 



it 



iR' v i d 

R^p f df 



R' 



(A.4) 



where the terms could be simplified by making use of 
Bessel's equation. 

We split the wave functions (A.2) into angular and 
radial parts, <j)j((p) and Xj(?)i respectively, and introduce 
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the following matrices of f- and ^-integrals 



Jo 



K 



r [a\ _ 
X jk - 



lb] 



jk 



iRip' 



dip (f> 



v=JM 



R, 



v 



f>j4>k, 

(A.5) 





T [a] 


T [b] 
x jk 


T [c] 


Vk—pj = 
Vk—fJ-j = ±m 
Mfc —Mi = ±2m 

otherwise 


£ m ™ 2 A m /2 
0(A m ) 
-m 2 A m /4 

-r- mm' \ \ 
H 2 ^fn s^ rn 

o(\ 3 m ) 




±rn A m /2 
O(A^) 


1 

A m 

A m /4 
1 \ \ 





TABLE III: Expansion of the integrals X [a \ l [b \ and l [c] for 
small perturbations. Non-zero values are only found if one of 
the conditions given in the first column is satisfied. 



Note that the do not depend on the deformation 
parameters (p, A) and that the 1^ only depend on A. 
Decomposing the matrices Ti and B in terms of these 
newly defined matrices we find 



Hjk - 



n 2 V [b] \l [a] - 



if 7, 



(A.6) 

Ik, 



B 



Ok 



3 k 

if 3 



k, 

p 2r R [ ^l otherwise. 



(A.7) 



These results can now be implemented to solve the gen- 
eralized eigenvalue problem numerically. Obviously, the 
larger the deviations from axial symmetry, the more ba- 
sis functions have to be used in order to achieve good 
convergence. 

Expansion for small deformations of a cylinder: 
Moreover, Eqs. (A. 6) and (A.7) also allow to study ana- 
lytically the effect of a small deformation of a cylindrical 
wire on the transverse eigenenergies and the lifting of 
the degeneracies through the breaking of the rotational 
symmetry. Hereby we can give a proof of the relations 
(21-23) which ensure that the stability matrix of a cylin- 
drical wire is always diagonal. 

Let the initial shape be deformed according to the ra- 
dius function (20). The deformation parameters are con- 
sidered to be small so that we can expand the integrals 

l up to second order in the A m . We find that they take 
non-zero values only if the quantum numbers pj and fj, k 
satisfy one of the following relations, 



Pk — Pj = in - m, 
\fi k - Hj\ = \m±fn\, 



(A.8) 
(A.9) 



(/i = 0) it is found that 



n 



(i) 

33 



E {0) B {1) 

3 33 



(A.10) 
(A.ll) 



EfB^- 



E 

k±3 



n 



jk 



E 



(o) 



E 



(o) 



where the matrix B in the equations above reflects the 
use of a non-orthogonal basis. From Eqs. (A. 6) and (A.7) 

and the first row of Table III we deduce H 
Bf) = 0, so that Ef ] 



(i) 



B {1) = 

33 



and 



P E i 



n 



(2) 
33 



* E 



n 



a). 

jk 



(o) 



(A.12) 



= E^ 



V 



2/to 2 



+ E 

iMfc — 1 = 



fV^fc -ti [c] 
2 'S'fe 



1 3 '^3k 



it 



Since E^ does not contain mixed terms oc \ m ^m (with 
m ^ m) we have hereby derived the relation (21). 

The case p > is slightly more complicated as it re- 
quires degenerate perturbation theory: In first order we 
have to solve 



K 



(1) /»« K?(l)\ 

31 jj I \ 3 1 33 / 



= o, 

(A.13) 



where 1 is the unit matrix, and i = (p,v), and j — 
(—p, v). Using the expansions given in Tab. Ill, we find 



where n is an integer and to and to are elements from 
the set IM of to- values appearing in the sum of Eq. (20). 
The results are summarized in Tab. III. 

We can now straightforwardly derive the energy levels 
Ej within perturbation theory of the generalized eigen- 
value problem. Therefore we expand all the quantities 
entering the problem up to second order in the perturba- 
tion and solve order by order. For non-degenerate states 



E 



i i = J ±(7i/p) 2 if 2/i = to, 
| otherwise. 



(Add) 



The energy change in second order degenerate pertur- 
bation theory is found to read 



E 



(2). 



f)A.15) 
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where we have used B\P — = 0, and the matrix is equal to the full result (A. 12) for the non-degenerate 

case. Therefore we conclude that relations (22) and (23) 
r V (n^-E^)(H^-E<^) hold. 

(2) 

Note that the quantity H u +Cu appearing in Eq. (A. 15) 
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